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Abstract 

We consider a system with N unit-service-rate queues in tandem, with exogenous arrivals of rate A 
at queue 1, under a back- pressure (Max Weight) algorithm: service at queue n is blocked unless its queue 
length is greater than that of next queue n + 1. The question addressed is how steady- state queues scale 
as N — > oo. We show that the answer depends on whether A is below or above the critical value 1/4: in 
the former case queues remain uniformly stochastically bounded, while otherwise they grow to infinity. 

The problem is essentially reduced to the behavior of the system with infinite number of queues in 
tandem, which is studied using tools from interacting particle systems theory. In particular, the criticality 
of load 1/4 is closely related to the fact that this is the maximum possible flux (flow rate) of a stationary 
totally asymmetric simple exclusion process. 
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1 Introduction 

In this paper we consider a system with N single-server queues in tandem, with exogenous Poisson arrival 
process at queue 1 and customers leaving after service in queue N. All service times have exponential 
distribution with unit mean. The waiting room in each queue is unlimited, but the system operates under 
the following back-pressure policy: service of queue n is blocked (stopped) unless its queue length is greater 
than that of "next- hop" queue n + 1. The system is stable as long as input rate A < 1. The main question 
we address is: Given load A < 1 is fixed, how steady-state queue lengths increase (scale) as N — > oo. 

Our main motivation comes from the fact that general back-pressure (BP) policies (sometimes also called 
MaxWeight), originally introduced by Tassiulas and Ephremides [12] and received much attention in the 
literature (cf . [121 [4] for recent reviews) , are very attractive for application in communication and service 
networks. This is due to the adaptive nature of BP policies - they can ensure maximum possible network 
throughput, without a priori knowledge of flow input rates. Key mechanism in BP policies, giving it the 
adaptivity (and the name), is that the "priority" of a traffic flow / at a given network node n is "proportional" 
to the difference between flow / queue lengths at node n and the next node on the route; in particular, 



unless this queue-differential is positive, flow / service at n can be blocked. This mechanism, however, has 
a drawback, namely it can lead to a large queue build-up along a flow route, since, roughly speaking, the 
queue "needs" to increase as we move back from a flow destination node to its source. Such "bad scaling" 
behavior of BP algorithm is emphasized in [3] , and is not very surprising (see also our Proposition [5] and 
bound (O). One approach to mitigate this scaling problem in practical systems is proposed in [2j[3], where 
it is suggested to "run" BP algorithm on virtual queues as opposed to physical ones. 

It is of interest to understand fundamental scaling properties of BP algorithm. In this paper, we address 
this problem for a simple system - a single flow served by iV queues in tandem. We show that, perhaps 
somewhat surprisingly, the scaling behavior of BP is not always bad. Namely, if load is below some critical 
level, A < 1/4 for our model, all queues remain uniformly stochastically bounded for all N. (In fact, we show 
that the stochastic bound has exponentially decaying tail; see Theorem 0Ji).) When A > 1/4, the queues 
increase to infinity with TV (see Theorem [TJii)). 

The problem of asymptotic behavior of queues as N — )• oo is essentially reduced (see Proposition [4]) to the 
behavior of the system with infinite number of queues in tandem. Such infinite-tandem model is within the 
framework of interacting particle systems [3 [8] - methods and results of the corresponding theory will be our 
main tools. As we will see, the criticality of load 1/4 is closely related to the fact that 1/4 is the maximum 
possible flux (average flow rate) of a stationary totally asymmetric simple exclusion process (TASEP). In the 
subcritical case, A < 1/4, a stationary regime exists such that only a finite (random) number of "left-most" 
queues can be greater than 1, the rest of the queues have at most one customer and the process "there" 
behaves like TASEP (see Theorems |8] and |9|) . In the supercritical case, A > 1/4, each queue grows without 
bound with time (see TheoremllOp. 

Infinite series of queues in tandem is a much studied model in the literature (cf. [llj and references therein) , 
under a variety of assumptions. In particular, |11) studies an infinite system under a class of blocking policies 
(which are different from the BP policy in our model), where blocking is caused by limited waiting space 
(finite buffer) in the queues; there, the phenomenon of critical load, below which the system is stable, also 
exists. 

The rest of the paper is organized as follows. Section [2] presents the formal model and main result. The 
"reduction" of our problem to the behavior of infinite system, and basic properties of the latter, are given 
in Section [3] The subcritical and supercritical load cases are treated in Sections |4] and respectively. In 
Section [5] we remark on more general input processes. 



2 Formal Model and Main Result 



Consider a series of N servers (sites), numbered 1, . . . , N, each with unlimited queueing room. A Poisson 
flow of customers, of rate A > 0, arrives at server 1, and each customer has to be served consecutively by 
the series of servers, from 1 to N; after service by iV-th server, a customer leaves the system. The service 
time of any customer at each server is exponentially distributed, with mean value 1; all service times are 
independent of each other and of the input process. To be specific, assume that each site serves customers 
in first-come-first-serve order - given Markov assumptions, this will not limit generality of results. 

Let us denote by Q„ (t), n = 1, . . . , N, t > 0, the queue length at n-th server at time t. The superscript N 
indicates the number of servers, which will be the parameter we vary. 

Consider the following back-pressure (BP) algorithm: site n is actually serving a customer (which is the 
head-of-the-line customer from its queue) at time t if and only if Qn(t) > Q^T^t). (We use convention 
that Q^l^t) = 0.) In other words, the service at site n is blocked unless the queue at site n + 1 is smaller. 

The random process = (Q ( n\t), n = l,2,...,N), t > 0), describing evolution of the queues, is 

a countable irreducible continuos-time Markov chain. Stability of this process - ergodicity of the Markov 
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chain - is guaranteed under condition A < 1 - this follows from well-known properties of BP algorithms (cf. 
[3]). Therefore the unique stationary distribution exists; we denote by Q^ N \oo) a random system state in 
the stationary regime. 

The question we address is how steady-state queues Qn (oo) grow (scale) as N —> oo, more specifically 
whether or not they remain stochastically bounded. It is easy to observe that if condition 

Q [ nHt) + 1 > maxQ ( k N) (t), Vn > 1, (1) 

holds for t = then it holds for all t > to as well; in particular, it does hold in the stationary regime. Thus, 
Q[ N) (oo) + 1 is an upper bound on all queues. Therefore, we can concentrate on the question of whether 
or not (oo) remains stochastically bounded as N — > oo. (If it does not, it is easy to see that, for any 
n, Q n N \oo) goes to infinity in probability.) We show that the answer depends on whether or not the input 
rate A is below or above the critical value 1/4. Namely, our main result is the following 

Theorem 1. (i) If A < 1/4, then there exist C\ > and Ci > such that, uniformly on N, 

P{Q[ N \oo)>r}<C ie - C2r . (2) 
(ii) If A > 1/4, then Q\ (oo) —¥ oo in probability as N — > oo. 

Statements (i) and (ii) will follow from Theorems[51and ll0l respectively, which concern with the corresponding 
infinite-tandem system. 

Note that the fact that tight uniform bound ([2]) cannot possibly hold for all A < 1 is fairly obvious. Observe 
that in the stationary regime the average rate at which customers move from site n to n+ 1, for any n < N, 
is A. Therefore, we have 

Proposition 2. If A < 1, then for any n = 1, . . . , N, 

P{QW(oo)>Q£0 (cx>)} = A. 

Then, using (QP, 

EQW(oc) - EQW (oo) > 2A - 1, n = 1, . . . , N. 

Thus, in the case A > 1/2, we have at least linear growth of the first queue expected value: 

EQW(oo) > {2X-1)N. (3) 



3 Basic Properties of Infinite- Tandem Queues under Back-pressure 
algorithm 

Consider a system just like the one in Section[2j except there is an infinite number of servers (sites) in tandem, 

indexed n = 1, 2, Arriving customers never leave - they just keep moving from site to site, to the "right" . 

We denote by Q n {t), n = 1,2, . . ., t > 0, the queue length at site n at time t, by Q{t) = (Q n (t), n = 1,2, . . .) 
the state of the entire system at time t. It will be convenient to assume that the phase space for each 
queue (site) n state Q n is the compact set Z + = Z + U {oo} with Z + being the set of non- negative integers 
and with metric, e.g., \e~ % — e~ J \. (We use convention oo — 1 = oo, so that the back-pressure algorithm is 
well-defined.) The state space of Markov process {Q{t), t > 0) is Z+, S = {1, 2, . . .}, with product topology; 
the process is formally defined within the framework of interacting particle systems (cf. Section 1.3 of [7], 
specifically Theorem 1.3.9). We will use some sligtly abusive notations: Q(-) for the process (Q(t), t > 0), 
and Q = {Q n > n — 1,2, . . .) S for elements of the phase space; and will denote ||Q|| = ^2 n Q n . 
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Throughout the paper we will also use the following representation of process Q(-), which is standard for 
Markov interacting-particle systems and is convenient for coupling (cf. [7J |S]) of different versions of the 
processes. First, assume that there is a site associated with each integer n G Z, not just positive n. (We will 
use this convention in Section 21) The underlying probability space is such that there is a unit rate Poisson 
processes II„ for each site n; these processes are independent from each other and of the input flow Poisson 
process. Then, if time r > is a point of the Poisson process II„, a customer moves from queue n to queue 
n+ 1 at t if Q n (r—) > Q n _|_i(r— ), otherwise the move is suppressed. 

The finite system of Section [3J with N sites, will be viewed as an infinite one, but with the modification that 
any customer reaching site N + 1 is immediately removed from the system, and with Q n (t) = for n > N 
(and n < 0). 

It is easy to check that an analog of ([T]) holds for the infinite system as well. Namely, if 

Qn{t) + 1 > sup Q k (t), Vn > 1, (4) 

holds for t = to > 0, then it holds for all t > to as well. 

Now we state basic monotonicity properties of the infinite system. (They are easily proved by contradiction, 
using coupling on the common probability space defined above.) The inequality Q < Q' is understood 
component-wise; the order relation Q ■< Q' means that both \\Q\\ and ||Q'|| are finite and 

k>n k>n 

Lemma 3. If Q(0) < Q'(0) [respectively, Q(0) ^ Q'(0)], then the processes Q(-) and Q'{ ) can be coupled 
so that Q(t) < Q'(t) [respectively, Q(t) < Q'(t)] for all t > 0. 

As a corollary of Lemma [3l we obtain the following 

Proposition 4. Consider the infinite system and the finite systems, for each N = 1,2,..., all with zero 
initial state (with all queues being 0). Then, 

(i) All corresponding processes can be coupled so that for all t > 0, 

Q«(t) < Q {2 \t) < . . .QW(t) < . . . < Q(t), 
Q (1) (t) d Q (2) W d ■■■Q {N) (t) d...dQ(t). 

(ii) Process Q(-) is stochastically non- decreasing in t (in the sense of < order), and process Q^ N \-) is 
stochastically non- decreasing in both t and N . 

(Hi) We have convergences in distribution 

Q^ N \t)^Q {N \oo), ast^oo, (5) 

Q (N) {oo) ^ Q(oo), asN->oo, (6) 
where the sequence in the left-hand side is stochastically non- decreasing and Q(oo) just denotes its limit, 

Q(t) => Q(oo), as t -> oo. (7) 

(iv) The distribution of the limit Q(oo) is a stationary distribution (namely, the lower invariant measure) of 
Markov process Q(-). 

(v) Condition Q) holds for all t > and for the stationary state Q(oo). 

Proposition 21 iii) of course implies that Q\ (oo) is stochastically non-decreasing, converging in distribution 
to Qi(oo) and therefore, to prove Theorem [Tl we can study the distribution of Qi(oo). (If it happens that 
<5i(oo) = oo, this implies Q^(oo) — > oo in probability as N — > oo.) 
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We will need one more monotonicity property, which is also a corollary of Lemma [3] Its meaning is very 
simple: if in addition to process Q(-) we consider another process Q'(-) which is constructed the same way 
as Q(-), but with some additional exclusions ("obstructions") on the movement of the customers, then Q'(-) 
will stay "behind" Q(-) in the sense of < order. 

Proposition 5. Consider a fixed realization of the process Q(-), which includes a fixed finite initial state 
Q(0) and realizations of the input process (at site 1) and of all processes II„, n > 1. Now, suppose further 
that in the realizations of processes ±1„, some of the points (jumps) are marked as "valid" (in an arbitrary 
way) and the remaining points are "invalid" . Consider another realization Q'(-), with the same initial state 
Q'{0) — Q(0), and constructed in the same way as Q{t), except invalid points of processes H n are "ignored" 
(cause no action). Then, 

Q'(t) r< Q(t), V£ > 0, 

and, in particular, Q'i(t) > Qi(t) for all t. 

4 Subcritical case: A < 1/4 

Suppose A < 1/4. We will construct a process Q'(-), coupled with Q(-) so that Proposition[5]holds path-wise, 
and such that we can obtain a stochastic upper bound on Q'i(t). 

Consider process Q(-), with zero initial state ||Q(0)|| = 0, constructed on the probability space described in 
Section [3] We will extend the probability space to define a stationary totally asymmetric simple exclusion 
process (TASEP, cf. Chapter VIII of [7]), with sites being integers n G Z, and particles moving to the "right" . 
Specifically, let us choose arbitrary (density) p < 1/2, such that p = p(l — p) > A. Let Y n (t) 6 {0, 1} denote 
the number of particles of TASEP at site n at time t > 0. We augment the probability space so that, 
independently of all other driving processes, at time each site neZ contains a particle, 5^(0) = 1, with 
probability p and does not contain one, Y n (0) — 0, with probability 1 — p. The movement of TASEP particles 
will be driven by the same Poisson processes H n , that drive process Q(-). (The exogenous input process at 
site 1 does not affect TASEP.) If time r is a point (jump) of II n associated with site n, then the particle 
located at n (if any) attempts to jump to site n + 1 - it actually does jump if site n + 1 is empty, and it stays 
at n otherwise. It is well known that if the initial state Y(0) has Benoulli distribution as defined above, 
then the TASEP process Y(-) is stationary (cf. Theorem VIII. 2.1 of [7]). The flux of this process, i.e. the 
average departure rate of particles from a given site, is p = p(l — p); the average speed of a given ( "tagged" ) 
particle is v — 1 — p. (For example, p = 1/2 gives the maximum possible flux p = 1/4 > A; this is where the 
condition A < 1/4 comes from: A needs to be less than the flux of a stationary TASEP.) 

The process Q'(-) has the same (zero) initial state as Q(-), and is constructed the same way as Q(-) except 
for an additional exclusion: a customer (particle) from queue (site) 1 cannot move to queue 2 at time r 
unless a particle of the TASEP jumps from site 1 to 2 at r. We now record basic properties of process Q'(-). 

Proposition 6. (i) Q'(t) < Q(t) and Q[(t) > Qi(£) for all t > (by Proposition^. 

(ii) For any n > 2 and any t > 0, Q' n (t) < Y n (t). In other words, at all sites to the right of 1, process Q'(-) 

stays "within TASEP"; in particular, there can be at most one particle in each site n > 2. 

(Hi) A particle jump from site 1 to 2 in the process Q'(-) happens at time r if and only if Qi(t—) > 1 and 

there is a jump of TASEP particle from 1 to 2 at time r. 

(iv) Process Q'(-) is stochastically non- decreasing with t. 

We know from Proposition [S] that Q[(t) is an upper bound of Q\(t). The behavior of queue length Q[(t) is 
such that it is initially zero, Q'i(0) = 0, the input process is Poisson with rate A, and the "service process" 
is the stationary process of TASEP particle jumps from site 1 to 2. We denote by A(tx, tz) and S(tx,t2) the 
number of points (jumps) of the arrival and service processes, respectively, in the interval (t\, £2]; WLOG we 
assume that these processes are defined for all real times, that is £i,£2 £ K, £1 < £2; their average rates are 
A and p = ES{t\, £2)/^ — £1), respectively. From the large deviations estimate given below in LemmaO it 
also follows that S{— s,0)/s — > p, s — > 00, with probability 1. 
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From classical Loynes constructions |10) it is known that the distribution of Q'i(t) is stochastically non- 
decreasing with t, and as t — » oo it weakly converges to the stationary distribution, which in turn is same as 
that of random variable 

Qi(oo) =sup[A(-s,0) -S(-s,0)]. (8) 

s>0 

Q'l(oo) is a proper random variable due to condition A < /i, which guarantees that the RHS of © is finite 
w.p.l. 

Thus, we see that, as t — >■ oo, (t) (and then Qi(t)) remains stochastically bounded by Q^(oo). Moreover, 
the large deviations estimates (Lemma [7J imply exponential bound on the tail of Q^(oo) distribution. We 
proceed with the details. 

The following fact is a known property of the stationary TASEP defined above. At time 0, let us consider 
the site with smallest index n$ > 2 that contains a particle, and tag this particle. (Obviously, uq — 2 has 
geometric distribution.) If we consider the point process of jumps of tagged particle in time interval [0, oo), 
it is a Poisson process of rate v — 1 — p (cf. Corollary VIII. 4. 9 of [7]). Therefore, the location H(t) of tagged 
particle at time t > is Hit) = 2 + Hi + Hzif), where Hi is geometric random variable with mean (1 — p)/p, 
H2(t) is Poisson r.v. with mean vt, and Hi and H2(t) are independent. From the stationarity of TASEP 
we also know that, at any time t, the total number G(n) of particles at sites 2, 3, n is simply the sum of 
n — 1 independent Bernoulli variables with mean p. Using "separate" large deviations estimates for H (t) 
and G{n), even though these two r.v. are not independent, we obtain the following 

Lemma 7. For any 5 > 0, there exist C3 > and C4 > such that 

¥{\S{0,t) -fit\ > St} < C 3 e~ C4t . (9) 

Proof. To prove bound 

P{S(0, t)<(jt- 8)t} < C 3 e~ C4t 
we can choose e > 0, small enough, so that 

P{S{0, t)<(p- S)t} < ¥{H(t) <(v- e)t} + P{G((v - e)t) <(p-e)(v- e)t}. 

Bound 

P{5(0, t)>(ji + S)t} < C 3 e~ C4t . 
is proved similarly. □ 

Theorem 8. Assume A < 1/4. There exist Ci > and C2 > such that 

¥{Q[{oo) >r} <Cie~ C2r ; (10) 

and then 

P{Qi(oo) >r}< Cie- C2r . 



Remark. Given the large deviations bound ((9|), the argument to prove (|10|) is quite standard (see [6j 
[5]). However, formally, [5] for example, requires a stronger condition, large deviations principle (LDP) for 
S(0,t)/t; we did not find this LDP result in the literature and it is not needed for our purposes. Hence, for 
completeness, we give a proof of the theorem. 

Proof. It follows from the definition ([8]) that for any fixed d > 
<2i(oo)< sup [A(-kd,0)-S(-(k-l)d,0)} = A(-d,0)+ sup [A(-kd, -d) - S(-(k - l)d, 0)] (11) 

fe=l,2,... fe=l,2,... 

Let us fix b > such that bX < 1, and for each r > we will choose d = br. Then we can write 
P{Ql(oo) > r} < P{A(-6r,0) > r} + ^ ¥{A(-kd, -d) - S(-(k - l)d, 0) > 0}. 

/c=2,3,... 
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If we fix S > small enough so that A + S < fi — S and (A + S)b < 1, we have 

P{Qi(°o) >r}< ¥{A{-br, 0) > (A + 5)br}+ 
\P{A(-kd, -d) > (A + 5)(k - l)d} + P{S(-{k - l)d, 0) < (p - 5){k - l)d}}. 

fc=2,3,... 

It remains to apply large deviations bounds on A(-) and S(-) (Lemma[7]). □ 
Let us define 

B(Q(t)) = mm{n =1,2... | Q n (t) = 0}, 

n 

B'(Q(t)) = min{n = 1,2... | ^ Q fc (t) < n}. 

fc=i 

B(Q(t)) is the "busy interval" of the process Q(t): to the right of, and including, site B(Q(t)) all sites have 
at most one customer and so the (instantaneous) evolution of the process follows the same "rules" as that of 
TASEP. The interpretation of B'(Q(t)) is as follows: starting from the state Q(t), we take a customer from 
the left-most site with 2 or more customers, and move it to the left-most empty site, and then repeat until 
all sites have at most 1 customer; then site B'(Q(t)) is the left-most empty site of the modified state. 

Obviously, B'{Q(t)) > B(Q(t)); and it is easy to check that Q'(t) < Q(t) implies B'(Q'(t)) > B'(Q(t)). 
Thus, B'(Q'(t)) > B(Q(t)). 

Theorem 9. Assume A < 1/4. There exist C5 > and Cq > such that, for all t > 

F{B'(Q'(t)) >n}< C 5 e~ c ° n , 

and then 

¥{B(Q(t)) >n}< C 5 e- Cen . 

Proof. The result is easily derived from the following facts: Q'i(t) is stochastically bounded by Q^oo); the 
distribution of Q[(oo) has exponentially decaying tail (Theorem [8]) ; process Q'(t) "stays within" stationary 
TASEP of density p at all sites n > 2; definition of B' . □ 

Theorem |9] illustrates in particular the fact that (when A < 1/4) the infinite-tandem system under BP 
algorithm in stationary regime is such that there only a "small" number of sites (from site 1 to site B(Q(t))) 
where queue can be greater than 1; while all sites to the right of B(Q(t)) have queue of at most one, and 
therefore the behavior of the process "to the right of B(Q(t))" is same as that of TASEP. 



5 Supercritical case: A > 1/4 

Note that if A > 1/2 we immediately see from ([3]) that EQi(oo) = limjv EQi (oo) = oo. Here we prove 
that, in fact, Qi(oo) is infinite w.p.l, under a weaker condition A > 1/4. The intuition behind our argument 
is as follows. Unless Qi(oo) = oo w.p.l, busy interval B{Q(t)) must be stochastically bounded w.p.l. Then, 
in stationary regime, all sites "far enough" to the right have at most one customer (particle) in them, and 
therefore the process "there" behaves as TASEP. However, a stationary TASEP cannot have flux greater 
than 1/4, while the flux of our process must be A > 1/4, a contradiction. 

Theorem 10. Assume A > 1/4. Then Qi(oo) = oo. (And then Q n (oo) = oo for all n > 1.) 
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Proof. The proof is by contradiction - assume Qi(oo) is finite with positive probability. It is easy to see that 
this is possible only if Qi(oo) < oo with probability 1. (Otherwise, using the facts that Qi(t) is stochastically 
increasing and weakly converges to Qi(oo), and using coupling, we could show that Qi(oo) = oo with prob. 
1.) This in turn implies that B(Q(oo)) (which is the weak limit of stochastically non-decreasing r.v. B(Q(t))) 
is also finite w.p.l. (Otherwise, we again could show that Qx(oo) = oo with prob. 1.) 

The stationary version of the process Q(-) (i.e. the one with stationary distribution equal to that of Q(oo)) 
we denote by Q(-). This process is such that: w.p.l condition (Q} holds for all t > 0, and therefore all Q n (t) 
for all t are uniformly stochastically upper bounded by Qi(oo) + 1 and then finite w.p.l; B(Q(t)) arc finite 
w.p.l (equally distributed) random variables for all t; the flux is equal to A, namely, ~RF n (t)/t = A for any 
t > and n > 1, where F n (t) is the number of customers arrivals at site n in interval (0, t]. Note that 
for n > 2, the arrival process F n (-) at site n is the departure process from site n — 1. This implies that 
increments of processes F n (-), n > 2, are stochastically upper bounded by the increments of independent 
Poisson processes II n _i(-). 

Consider space-shifted processes {[T m Q](-),[T m F](-),[T m Il](-)}, m = 1,2,..., where [T m Q]i(t) = Q m+i (t), 
i = 1,2,..., t > 0, and [T m F]i(t) and [T m II]j(t) defined similarly. For each m this process is such that 
[T m <3](')] an d the increments of [T m F}(-) and [T m IT](-) are stationary. (Note that process [T m Q](-) is not 
Markov.) This process is still well-defined if we assume that each component [T m Q]i(t), [T m F]i(t) and 
[T m Il]i(t) takes values in (non-compact) space Z + with the usual topology (because we know that they are 
finite w.p.l), and with corresponding product topology in the space of values of the entire process. Note 
that, since B(Q(t)) is finite w.p.l and condition Q holds, we have 

lim P{sup[T m Q]„(i) < 1} = 1, Vt > 0. (12) 

m— yoo n>\ 

Then, using properties of process Q(-) described earlier, in particular the fact that increments of F n (-) are 
bounded by those of n„_i(-), it is easy to see that a process consisting of any finite subset of components 
[T m Q]i(-), [T m F]i(-) and [T m II]j(-) is tight (cf. Theorem 15.6 in [1]). Consequently, there exists a subsequence 
of {m} along which the shifted process converges in distribution to a process {Q(-), F (■),!!(■)}, which has 
the following (easily verifiable) structure and properties: 

(a) Q(-) is stationary, with flux equal A; 

(b) Q n (t) < 1 for all n and all t; 

(c) the movement of customers between sites is driven by independent, unit rate Poisson processes IT(-), 
according to BP algorithm rules; 

(d) By (b) and (c), Q(-) is a TASEP (with the "exogenous" arrivals at site 1 forming a stationary process, 
not independent of the "rest of the process"). 

Consider the following projection of process Q(-). All particles arriving at site 1 after time and the particle 
located at site 1 at time (if any), we will call "new" particles, while all particles initially present at sites 
n > 2 are "old". Let Q*(-) denote the process "keeping track" of new particles in Q(-), namely Q£(i) = 1 if 
there is a new particle located at site n at time t, and Qn(t) = otherwise. The flux of process Q*(-) from 
site 1 to site 2 is obviously equal to the flux of Q(-), which is A. We will compare Q*(-) to the following 
TASEP Q**(-), coupled to it - with the same Poisson processes driving movement between sites. The initial 
state of Q**(-) is: Qf(0) = 1 and Q**(0) = for n > 2. By definition, Qf*(t) = 1, that is if at any time a 
particle moves from site 1 to 2, it is immediately replaced at site 1 by another particle. Using the path-wise 
monotonicity properties given in Section |3l it is easy to see that Q*(t) ^ Q**(t), which implies 

FS(t) < F?(t), t>0, 

where F%(f) and i ? 2 **(^) are the numbers of particle arrivals in (0,t] at site 2 in the processes Q*(-) and 
Q**(-)) respectively. Then, 

hminf5£m> lim ^ = A . ( 13 ) 

t— >oc t t— ^oo t 

The TASEP Q**(-) is a special case of one of the processes studied in [5]. It is known (see Theorem 1.8(a) 
and Theorem 1.7(b) of [5]) that the distribution of Q**(t) converges to a stationary distribution, with 
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the corresponding stationary process having flux 1/4. This means that lim^oo KF£*(t)/t = 1/4, which 
contradicts (|13p since A > 1/4. Proof is complete. □ 



6 Remark on more general input processes 

The Poisson assumption on the input process is adopted to simplify exposition. We belive our main results 
can be easily generalized for the case of a stationary ergodic input process A(-), as long as large deviations 
bound 

¥{\A(fi,t) - At | > St} < C z eT Cl \ (14) 

analogous to (j9]), holds for any 5 > 0. Moreover, if (fT4|) does not hold, and we only have the ergodicity of 
A(-), the uniform stochastic boundedness results of Theorem [5] (and then Theorem HJi)) and Theorem |H] will 
still hold, except the bounds are proper (finite w.p.l) random variables, not necessarily with exponential 
tails. 

Acknowledgement. I would like to thank Yuliy Baryshnikov for extremely helpful discussions throughout 
the course of this work. 
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